Berry's phase for large spins in external fields 
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It is shown that even for large spins J the fundamental 
difference between integer and half-integer spins persists. In 
a quasi-classical description this difference enters via Berry's 
connection. This general phenomenon is derived and illus- 
trated for large spins confined to a plane by crystalline elec- 
tric fields. Physical realizations are rare-earth Nickel Boro- 
carbides. Magnetic moments for half-integer spin (Dy"^"'', 
J = 15/2) and magnetic susceptibilities for integer spin 
(Ho3+, J = 8) are calculated. Experiments are proposed to 
furnish evidence for the predicted fundamental difference. 
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The common belief that large spins are equivalent to 
classical charged gyroscopes was strikingly undermined 
by HaldaneEl who indicated that the groundstate and 
low-lying states of magnetic chains consisting of local- 
ized spins J are different for integer and half-integer J 
even if J is large. Here we consider a similar effect for 
individual large spins placed in an external fields such as 
the crystalline electric field (CEF) and the external mag- 
netic field. Its origin can be traced back to a geometric 
phase occurring on passing from quantum mechanical to 
quasi-classical spins. 

The general Hamitonian of a localized spin associated 
with a magnetic moment can be written as a function of 
its components: 



Hs - /(J) - hJ 



(1) 



where f{J) is an arbitrary function of J satisfying only 
two requirements: it is a Hermitian operator and it is an 
even function of J: f{J) = f{—J). The latter require- 
ment is equivalent to time-reversal symmetryB. In this 
paper we will focus on localized moments which are es- 
sentially confined to a plane. The confinement is due to 
the effects of the crystalline electric field. The important 
degree of freedom is rotation about the normal vector of 
a certain plane which wc choose as z-direction. It is then 
appropriate to introduce the azimuthal angle (f and its 
conjugate momentum Jz as canonical variables. Setting 
fi, = 1 we may use 



J. 



.d_ 
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Jx = \/ J{J + 1) — J| cos ip 
Jy ~ \/ J(J + 1) — J| sin if . 



(2a) 

(2b) 
(2c) 



Eqs. (||) generally have a symbolical meaning since the 
operators on the right hand sides should be ordered to 
guarantee that J+ = Jx] Jy — Jy ^-nd that the canoni- 
cal permutation relations are satisfied. However, at large 
J the non-commutativity is small and we do not need 
to bother with the symmetrization. Inserting (^ into 
the Hamitonian (|l|), the problem is reduced to the solu- 
tion of a Schrodinger wave equation with, in general, a 
complicated differential operator. 

So far, no difference between integer and half- integer 
spins occured. The difference resides in a different global 
phase behavior. For large spins J one may pass to a 
quasi-classical description via the overcomplete set of co- 
herent states -f/jfj with 



JlV'n) = Jn 



(3) 



Moving the spin induces a certain orbit of the tips of 
the unit vectors n on the unit sphere S"^ . If this orbit is 
closed there is no difference in the purely classical picture 
between starting point and end point. Due to Berry's 
connections, however, this is not the whole story. In 
order to be quantitative wc first have to fix the phases of 
the ijjfi in (H) . |-,The natural choice is to fix the phase of 
ijjs and to takeu 



ijjft = e-iip{iJz(p) e-npiiJ^O) exp(-i J^.ys)^'? 



(4) 



where (if, 0) are the angles characterizing the unit vector 
ft (i* being the unit vector in z-direction). Berry's con- 
nection is given for the single, non-degenerate state V'a 
(abehan caseSlafl) hy A = {ipfi\iV\ipfi). Since n is con- 
fined to the unit sphere, spherical coordinates are most 
appropriate and only the Ag and the A^ components 
matter. One finds 







(5a) 



Ag = {^Pn\^^/{^0M-t) = -^J^m 

A^ = {^Pn\^sin{ey'^/{^^)\^bn) 

= sin(0)-i(V'?| J.(l - cos(0)) + sm{9)Jy\i^^) 

= J(l-cos(6'))/sin(6i) . (5b) 

This connection then gives rise to the geometrical phase 
exp (i AdP) along the path 7. Since we are interested in 
this work in motion in the xy-plane only, we have 9 = Tr/2 
and thus the phase is exp(i JAi^s). As was to be expected 
from the physical origin of this phase it can be put to zero 
locally by an appropriate gauge. Globally, i.e. for com- 
plete tours of A(p — 2tt, one sees that no effect occurs for 
integer J but that for half-integer J a factor -1 applies 
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which cannot be gauged away One way to account for 
the phase behavior is to use ( Pa|) and aniiperiodic bound- 
ary condition for half-integer spins and periodic bound- 
ary condition for integer spins, respectively. A more ele- 
gant way is to stick to the connection A. This means we 
use 



J. 



(6) 



instead of (Pa|). Since, however, a change of gauge 
can alter A^p by any integer value one may use ( |2a| ) 
for integer spins. For half-integer spins we use Jz — 

-i{d/idip)-i/2). 

Thus, despite of the large value of J, the low-lying 
states of the Hamiltonian (|l]) for integer and half- 
integer spins are fundamentally different. In the ab- 
sence of an external magnetic field all stationary states 
of a half-integer spin are doubly degenerate (Kramers 
degenerac yl3). This means that, in analogy to the linear 
Stark effect, the groundstate may be characterized by a 
finite magnetic moment. To the contrary, the ground- 
state of an integer spin is non-degenerate for sufficiently 
low crystalline symmetry. Therefore, the magnetization 
in the groundstate is exactly zero and depends linearly 
on the magnetic field as long as the field is weak enough. 

To be more specific, we consider an important appli- 
cation of these general ideas to the ions of the rare earth 
elements Ho^+ and Dy^+. The triply charged ions are 
decisive for the magnetic moments in the compounds 
i?Ni2B2C (R stands for a rare earth element) the prop- 
erties_Q£_which attracted much attention in the last few 
yeargOilJ. The ion Ho^+ has 10 electrons in the 4f-shell. 
According to Hund's rule, it has the total spin 5* = 2, 
the orbital moment L — 6 and the total moment J — 8. 
The corresponding numbers for Dy'^"'" (9 electrons in the 
4f-shen) are: 5'=|,L = 5, J=^. Thus, both values 
of J are rather large and close to each other. 

All the compounds (it! =Y, La, Ho, Dy, Tm, Tb, Er, 
Yb) crystallize as perovskites with|-|the rare earth ions 
forming a tetragonal centered latticeO. The magnetic mo- 
ments of the Ho and Dy compounds are confined presum- 
ably in the a6-plane thus realizing the situation discussed 
above. The simplest crystal electric field (CEF) Hamito- 
nian Hqef displaying tetragonal symmetry reads: 



CEF 



2 ' 



2b{Jt + Jt). 



(7) 



with a, & > 0. Quartic (and higher) terms in Jz are ne- 
glected since we assume Jz ^ Jx,Jy For the same reason 
one can introduce the coordinate in a slightly simpli- 
fied way (compared to (^): Jx = J cos if, Jy = J simp, 
and Jz = —id/{dip) + A^ with A^ = for integer J and 
A^ = 1/2 for half- integer J. Together with the magnetic 
field contribution, the total Hamiltonian becomes: 



H 



iA. 



a f d_ 
2 

-hJ cos{ip — iph) 



-J''(3 + cos(4v3)) 



(8) 



where iph is the angle the magnetic field forms with the x- 
axis. The effective field h is the magnetic field times iWb- 
The Hamiltonian (||) was found in a previous workllHl in 
which it was also analyzed for integer J. 

It will be shown later that the CEF constants a and 
bJ^ are of the same order of magnitude. Due to the large 
value of J the potential energy has very deep minima near 
the points ip = (pi = ln/2,l e {0,1,2,3}. Let us denote 
the oscillatory states localized near each value of as 
1^). To be more specific, each |^) is the groundstatellj of 
Hi = -{a/2)d^/{dip^) + Ui{ip) with 



+ cosi-iip)) 

2bJ* 



for \(p~m < 
otherwise 



(9) 



Without loss of generality we assume that the corre- 
sponding eigenenergy is zero. Neglecting the overlap be- 
tween different |^), we find that the energy level is four- 
fold degenerate. The overlap lifts partly the degeneracy 
even in the absence of magnetic field. In complete anal- 
ogy to-the derivation of the dispersion of a tight-binding 
modellj we define the hopping matrix element 



w 



-(l\H - Hi+i\l + l) > . 



The hopping part of effective Hamiltonian is 

Hu, = -wiC{a) + C+ia)) 



(10) 



(11) 



where C is the unitary rotation operator which induces 
\l) — > 1^ -I- 1) exp(iQ;) and |3) |0)exp(za). For integer 
spin we have a = 0, for half-integer spin we use a = 7r/4 
resulting from 



/•(i + l)7r/2 

exp(iQ!) ~ exp I i / A^pdip 

^ JItt/2 



(12) 



This is the direct effect of the connection A^ in absolute 
analogy to the Peierls phase in tight-binding models in 
magnetic fields. For what follows it is important to note 
that w is exponentially small if the wells at ^pi are deep 
enough such that the ground states |Z) are well, i.e. ex- 
ponentially, localized. To estimate w we use the ansatz 
\l) cx exp(-| /;^^2 \/2(?7i(</3) - C/min)/adv3|) which is mo- 
tivated by the groundstate for the harmonic potential 
close to the minima and its natural extension in a WKB- 
type approach. The main contribution to the (^-integral 
for w comes from the vicinity oi ip — -k /A for I — and 
leads to 



w oc exp 




(13) 



The eigenstates and eigenvalues of ( |llD are easily read 
off since we deal with a translationally invariant, c? = 1, 
four-site tight-binding model. Thus the eigenstates are 
characterized by some momentum k G {0, ±7r/2, tt} 



1 ^ 



(14) 
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The corresponding eigenenergies are 
Ek — ~2w cos(fc + a) . 



(15) 



Thus for integer J (a = 0) there is a non-degenerate 
groundstate (fc = 0,Ek — —^w), two degenerate excited 
states (fc = ±7r/2, Ek = 0), and the highest excited state 
(fc = TT, Ek = 2w). For half-integer J (a = 7r/4) both the 
groundstate and the excited state are doubly degenerate 
with (fc = Q,-TT/2,Ek = -V2w) and (fc = TT/2,n,Ek = 
V2w). 

The main difference between integer and half-integer 
spin resides here in a different degeneracy of the ground 
states. Physically this difference becomes manifest, for 
instance, when a magnetic field h is applied. We first 
consider an in-plane field as in (^). For well- localized 
states \l) such a magnetic field is site-diagonal with ma- 
trix elements 



Hh = /iJcos(/7r/4 — (fih) 



(16) 



The eigenvalues of H^+Hh can be given analytically. For 
integer spins and half-integer spins one finds, respectively 



E, 



hint 



± 



± 



\ 



2w2 



/l2 



± 



2w2 
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hfhl (17a) 



2^2 



±h\ 2w'^ 



T 



h? ■ 



(17b) 



where we used /i as a shorthand for hJ . From (y_7|) one 
obtains for the groundstate energies in the limit of small 
magnetic fields 



Eint 
-Ehint 



~2w - / {Aw) 
-V2w - hJ/2. . 



(18a) 
(18b) 



As expected the correction is quadratic in the non- 
degenerate, integer J case, but linear in the degenerate, 
half-integer J case (as in the linear Stark-effect). So we 
have for integer J a finite T = susceptibility 

Xpianc = -g'fiUdE)y{dh^) - .gVy V(2^) • (19) 

For half-integer J the system shows a finite magnetic 
momentum 



-gfiBdE/{dh) = .g/ie J/2 



(20) 



leading to a Curie susceptibility diverging for T 0. 
This constitutes an essential difference between integer 
and half-integer spin. For large magnetic fields, how- 
ever, the difference vanishes since w becomes unimpor- 
tant. Asymptotically one obtains for hJ ^ w 



E fa ±hj cos{iph)', ±hj sin{iph) 



(21) 



irrespective of whether J is integer or half-integer. This 
equation implies that the corresponding magnetic mo- 
ments at saturation are directed along one of the four 
easy axes. This intermediate asymptotics is valid in the 
range w ^ hJ <C 6 where our tight-binding treatment 
stays valid. At significantly larger h the saturation mag- 
netic moment is directed parallel to the magnetic field. 
Comparing (|2^) and (|l]) we, find that the finite mag- 
netic moment at low magnetic field for half-integer J is 
precisely half of its saturation value (9.8/xb for Dy'^+; cf. 
10/iB for Ho'^+ with integer J). 

Let us consider now the action of magnetic field along 
the z-axis. To do so we rewrite 



a J 2a 



id 

dip 



A, 



2a 



(22) 



^From this we infer that the z-axis magnetic field adds a 
constant to the Hamiltonian and acts as if the connection 
is changed in the manner A'^ = A^p — h/a. Insert- 
ing A'^p in (|l^) shows that the Peierls phase a in ( pT| ) is 
changed like a ^ a — 7r/i/(2a). The effect of the Peierls 
phase change on the eigenenergies and on the groundstate 
energy in particular is easily found in (pa). So, for integer 
J, we have E = — 2it;cos(7r/i/(2a)) —~h?/{2a) and from 
this Xaxis — — 'k'^w / {20^). Due to the exponential 
smallness of w (see ([T^)) the second term is negligible 
compared to the first one and we find Xaxis — a^^ for 
integer J. Note that Xaxis is exponentially small com- 
pared to Xpianc in ( p^ ) which proves the consistency of 
our treatment for which we assumed that the moments 
are essentially confined to the plane. For half-integer spin 
we consider E = —2wcoa{Trh/{2a) — 7r/4) — h? /{2a) and 
derive /i = —gfiBdE/{dh) = g^Qw/{\/2a) for the local 
moment in z-direction. We obtain again a moment in 
z-direction which is exponentially small compared to the 
one in the plane ( pO| ) consistent with the outset of our 
theory. We observe that the difference between integer 
and half-integer J is also visible in the magnetic proper- 
ties perpendicular to the easy plane. 

The quasi-classical treatment implying the splitting of 
each oscillatory level into a quadruplet is valid provided 
that the level spacing between the localized oscillatory 
levels io is much larger than w. It is certainly incorrect 
for energies E close to the maximum potential energy 
C^max = 26 J2. One should have a rough idea on how 
many quadruplets our treatment can be expected to ap- 
ply to. In the cases of interest J = 15/2, 8 the total num- 
ber of states is 16 or 17. Thus one can expect that one or 
two quadruplets are well described by the effective four- 
site tight-binding model. This conclusion is confirmed 
by direct diagonalization of the 17x17 matrix in a model 
crystal fieldcJ. In order to have an estimate for the ratio 
b/a we compare the total number of states 2 J -f 1 with 
the number of levels N in the four wells of the poten- 
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tial U{(p) = —{bJ'^/2){3 + cos(4i^)) in the quasi-classical 
approach 



N 



7r/4 



v/2m(L/inax - U{ip))dip 



326J4 



(23) 



The area in phase space of the maximum classical orbit 
is divided by 27r?i to obtain an estimate for the number of 
states. Here fi' jm is set to a. Equating N to 2J-f 1 k, 2 J, 



we find bJ /a 



1.23, i.e. they are of the same 



order of magnitude. 

Let us discuss the experimental consequences of the 
difference between integer and half-integer J described 
above. The clearest manifestations would occur in dilute 
alloys of Ho and Dy in the compounds Lui_2:Ho2;Ni2B2C, 
Lui_2;Dy2;Ni2B2C or in similar alloys with Y in place of 
Lu. Both Lu'^^ and Y'^+ have zero magnetic moment. 
The direct check of our prediction would be the measure- 
ment of EPR spectra for these alloys at low temnerature 
T < w (at about 2K according to the estimate inll£l). We 
expect that the lowest four levels are split as described 
above into singlet, doublet and singlet for Ho. The res- 
onance frequency will not depend much on the magnetic 
field as long as /i J <C w. In the case of Dy-alloys, how- 
ever, the EPR resonance frequency is proportional to the 
magnetic field as long as this is small due to the finite 
magnetic moment. O ne can als o check the nonlinear de- 
pendence of energies ( |l7a| , |l7b| ) on the magnetic field. 

Another option is to measure the magnetic field at a 
nucleus by NMR. The magnetic field at a nucleus must be 
essentially zero for Ho-alloys and non-zero for Dy-alloys. 

Finally, one can also measure the spin magnetic sus- 
ceptibility of the alloys at low temperatures. It should be 
constant (^/ib J)^/(2w) per ion in the Ho-alloys. In the 
Dy-alloys, an in-plane Curie susceptibility {giJ,)^/{3T) « 
(25(/^/x|)/(3T) per ion is expected due to the magnetic 
moments. Unfortunately, the paramagnetic effect will 
be strongly masked by the Meissner susceptibility xm = 
— 1/ (47r) due to superconducting currents. One can, how- 
ever, observe a remarkable weakening of the Meissner sus- 
ceptibility at the temperature ^ O.IK for concentrations 
of Dy of about 1%. At slightly lower temperatures the 
Dy-alloys must go over into a spin-glass state. No such 
transition can be observed in the Ho-alloys. 

In conclusion, we have shown that the well-known dif- 
ference between the magnetic properties of ions of in- 
teger or half-integer total magnetic moments J due to 
the absence or presence of Kramers degeneracy is cap- 
tured in the quasi-classical limit of large J by a Berry 
phase. In the tetragonal environment studied here, in- 
teger spins display no finite magnetic moment in the 
groundstate, whereas the half-integer spins display a fi- 
nite magnetic moment of half their saturation value. The 
difference is due to Kramers degeneracyl3 necessarily as- 
sociated with half-integer spins. Any quasi-classical ap- 
proach has to take the geometric Berry phase into ac- 
count in order to capture the essential difference between 
integer and half-integer spin. The difference should be 



manifest in the spectrum of EPR, NMR and in the mag- 
netic susceptiblity of dilute the alloys Lui_2;Ho2,Ni2B2C 
and Lui_2;Dya;Ni2B2C. 
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